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Interference of single-photon wave packets at a beam splitter usu- 
ally leads to an anticorrelation of the light intensity in the two 
output ports of the beam splitter. The effect may be regarded as 
"bunching" of the photons at the beam splitter and has widely been 
interpreted as a result of quantum mechanical interference between 
the probability amplitudes of indistinguishable bosonic particles. 
Here we show that when the wave packets are sufficiently distorted, 
then the opposite behaviour is observed, i.e., simultaneous clicks of 
the photodetectors in the two output ports are favoured, which 
may be regarded as "antibunching" of the photons at the beam 
splitter. 
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1 Introduction 



The study of the interference behaviour of quantized light has not only been 
of fundamental interest, but has also been important for applications, such as 
the use of nonclassical light in highly sensitive interferometry. In particular, 
two-photon interference experiments have offered new possibilities of study- 
ing fourth-order interference phenomena of single-photon wave packets [1—5]. 
An interesting phenomenon is the anticorrelation of the photoelectric counts 
recorded in the two output channels of a beam splitter [6] or an optical fibre 
multiport [7]. The dependence of the coincidence events on the optical delay 
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between the two interfering wave packets can be used to determine the photon 
tunneling time through an optical barrier [8-10]. 

In order to explain the anticorrelation (bunching) behaviour of interfering 
single-photon wave packets in the output channels of a beam splitter, it has 
been argued that when two photons arrive at the beam splitter simultane- 
ously, then the outgoing state is obtained by superposition of probability am- 
plitudes, because of the indistinguishability of the photons. Such an interpre- 
tation could imply the conclusion that the interference phenomenon observed 
is a pure quantum effect. On the other hand it has been shown that the ef- 
fect also appears for independent classical light fields [11]. It may therefore be 
expected that classical and nonclassical light fields give rise to quantitatively 
rather than qualitatively different interference structures of the coincidence- 
event statistics. In this context, the question arises of what is the effect on the 
interference fringes of correlated incoming beams and whether or not the in- 
terference fringes enable one to distinguish between classical and nonclassical 
light. 

In this paper we show that when single-photon wave packets whose superposi- 
tion gives rise to anticorrelations in the coincidence-event statistics in the two 
output ports of a beam splitter are sufficiently strong distorted by a multilayer 
dielectric slab and a superposition of the distorted wave packets is measured, 
then correlations can be observed in place of anticorrelations. Such a behaviour 
may be regarded as a kind of "antibunching" , which would be expected for 
fermionic particles rather than photons. Clearly, identifying the distorted wave 
packets with single photons can be quite questionable, as has been stressed 
recently [10]. Moreover, correlations can only be observed when the two in- 
coming wave packets are produced by the same source and the overall field is 
in a correlated state. We analyse possible interference structures and discuss 
criteria for deciding from the observed interference fringes whether the incom- 
ing fields are classical or not and whether they are correlated or independent 
of each other. 

The paper is organized as follows. In Sec. 2 the underlying formalism for 
describing the propagation of quantized light through multilayer dielectric 
structures is outlined. In Sec. 3 a Mach-Zehnder interferometer is considered 
and the interference fringes of the output coincidence events are determined in 
dependence on the distortion of the beams in the two arms. Finally, a summary 
is given in Sec. 4. 
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2 Field quantization and input-output relations 



To describe the interaction of quantized radiation with (linearly responding) 
dielectric matter, we use the quantization scheme developed in [12-14], which 
is based on the determination of the Green function of the classical propagation 
problem. In particular, we consider linearly polarized light propagating in x 
direction through a dielectric plate consisting of N = M — 2 layers (M > 3). 
The operator of the vector potential can then be given by [14] 
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1 if Xj-i < x < Xj and Xj(x) 



otherwise, ej{uj) — n^{uj) with 
Tij(u) =Pj(u) +ijj(u) is the complex permittivity of the jth layer, and A is the 
normalization area in the yz plane [j = 1 and j — M refer to the left- and right- 
hand side free-space surroundings of the plate, respectively, with 6i(u) = 6m(w) 
= 1]. The quasi- mode operators a,j±(x,u) inside the plate are associated with 
(damped) waves propagating to the right and left, respectively. Outside the 
plate they are independent of x and obey the familiar bosonic commutation 
relations for photons in free space, 
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where the output operators ai_(a;), a M+ {uj) are related to the input operators 
ai+(u>), clm-{w) an d bosonic operators of the plate excitations, g±(co), as 





9±(w),9t(w'j\ =5(u-u'). 



\ 



(5) 
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In Eq. (5), T and A, respectively, are the characteristic transformation and 
absorption matrices of the multilayer plate [14] . Note, that the elements of T 
and A are not independent of each other, but they are related to each other 
such that the consistency of the quantization scheme is ensured. 
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3 Fourth-order interference 



Let us consider an experimental scheme of the type studied in [8] (Fig. 1). 
Pairs of correlated (single-photon) wave packets produced in a parametric 
down conversion process are combined by a lossless 50%:50% beam splitter 
(BS) and the time-integrated output coincidences 

R(s)=eJ Ah J dt 2 (E^\t 1 )^-\t 2 )E^\t 1 )E^\t 2 )) (7) 

are measured (the detectors PDi and PD2 in Fig. 1 are assumed to have equal 
efficiency £). In contrast to [8], here each of the two wave packets is allowed 
to pass a dielectric barrier (DB 1 and DB 2 in Fig. 1) before impinging on the 
beam splitter. The position of a prism (P) in one arm of the interferometer 
is adjusted such that the additional path (2s) compensates for the difference 
between the propagations of the two wave packets trough different barriers. 

To calculate R(s), Eq. (7), we first relate the fields -E ± (ti) and -E ± (t 2 ) at the 
two detectors to the fields produced by the source, applying the input-output 
relations (5) and assuming that the barriers are in the ground state (zero- 
temperature limit). In a quantum- mechanical description we then perform 
the averaging (• • •) in Eq. (7) assuming that the down-conversion photons are 
initially prepared in an entangled state 

00 f! 

|*) = J dtoa(to) Jdwf(w)f(n-u)al(w)ah(to-u)\0), (8) 


where the subscripts I and II refer to the two arms of the interferometer (see 
Fig. 1), and a(Q) and f(u>), respectively, denote the band-shape functions of 
the laser and the down- conversion photons, f(ou) being centered at Q/2. After 
some calculation we derive 

R(s) = 27r 2 A/T [K - Ms)} , (9) 



Af= ^h/(4irceA), where 

00 Q 
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is independent of the position of the prism, s, and the position-dependent part 
Ki(s) reads as 



Tg>(u)f(u)T£>(n-u)f(n-u) ' (10) 
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K!(s) = jdna 2 (n) Jduouo{tt-uo) \f{uo)f{tt-uo)\ 2 



x e^-^TS\u)TS(Q - c)T 2 f*(Q - u)lg>(u). (11) 

In Eqs. (10) and (11), T^ lY) = 1,2) are the elements of the characteristic 
transformation matrices T^ 1 / 11 ) of the two dielectric barriers (for analytical 
results for lossless and lossy Bragg mirrors, respectively, see [15] and [14]). 

Let us compare this result with that obtained when a classical light source is 
considered which in a complete random way produces correlated pairs of wave 
packets. Straightforward calculation yields [16] 



R(s)=4n 2 e [G -G 1 (s)}, 
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Here, (• • -) c i denotes a classical averaging, the complex-valued Fourier com- 
ponents of the emitted (down-conversion) fields in the two arms of the inter- 
ferometer, E®(u) and E^ n \uj), being the random variables. Note that from 
Eqs. (13) and (14) the relation 



IGiOOl < i 
Go ~ 2 



(15) 



can be derived [16]. 

To perform the classical averaging, the (joint) probability distributions for 
E®(u) and E^(u) must be specified. To compare the quantum regime as 
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given by the state in Eq. (8) with a classical regime that is as similar as 
possible to the quantum regime, we assume that the probability distributions 
are such that the frequencies and phases of the laser fields and the down- 
conversion fields are related to each other as 8(0, — oo — uo') and 5(y? L — </?i — ipu) 
(</?l is the laser phase and ipi, ipu are the phases of the two down-conversion 
fields). In the spirit of close correspondence between classical and quantum 
regime we then may assume that Eqs. (13) and (14) can be rewritten as 



oo 

G = 2j duu(Q-u) \T^{u)f{u)T^\Q-u)f{Q-u) 



OO 

+ J duu(0-u) \T${u)f(u)T§(n-u)f(n-u) 





oo 



+ J dujuj{Vt-uj) \T$P(u)f(u)T$P(n-u)f(£l-u)\ 2 (16) 



and 
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where f(uj) is the same band-shape function as in Eq. (8). 



(17) 



The classical relations (16) and (17) and the quantum-mechanical relations 
(10) and (11) are seen to mainly differ in the second and third terms in 
Eq. (16). Physically, these terms are related to coincidence events caused ei- 
ther by the field in the arm I or the field in the arm II, since a classical field 
in one arm of the interferometer can always be decomposed into two parts 
by means of the beam splitter. The corresponding coincidence events are of 
course independent of the prism position s and hence the fringe visibility 



V 



(18) 



diminishes. Clearly, in the quantum regime a single photon in one of the arms 
cannot give rise to a coincidence event. 

Let us first assume that the barriers in the arms of the interferometer are 
absent and the prism is translated such that the two wave packets do not 
overlap at the beam splitter, i.e., s—^oo. As expected, there is no s-dependent 
interference effect, because the interference terms in Eqs. (9) and (12) vanish, 
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i.e., Gi(s) — > and Ki(s) — > 0, as it can easily be seen from Eqs. (11) and 
(17), respectively. In the quantum regime single-photon wave packets impinge 
separately on the beam splitter, where they are either reflected or transmitted 
with 50% probability in each case. In the classical regime the incident wave 
packets correspond, in a sense, to multiphoton states, because there is no 
classical analogue of a single-photon state. Hence, a wave packet that impinges 
on the beam splitter is split into two parts of equal intensity, which implies 
that coincidence events are preferably observed [see the second and third term 
in Eq. (16) - terms that are not observed in Eq. (10)]. 

To study the interference behaviour, let us introduce the normalized time- 
integrated coincidence events R^(s) = R(s)/R(oo) = [K — ifi(s)] / K [or = 
[Go-Gi(s)] I Gq in the classical regime]. From Fig. 2 it can be seen that in the 
quantum regime (s) decreases with decreasing value of s (compared to s — > 
oo) and attains a zero-value minimum. This behaviour can be explained from 
the argument of the indistinguishability of the two (now interfering) photons 
- there is no "which- way" -information available at the detectors, because of 
the wave combination at the beam splitter. The minimum is observed for 
equal optical lengths of the two paths I and II through the interferometer 
(i.e., s = in the figure). It should be pointed out that a qualitatively similar 
behaviour is also observed in the classical regime, which can be explained using 
standard arguments of classical optics. The difference between the quantum 
and classical regime is that in the latter the fringe visibility is substantially 
reduced. Recalling Eq. (15), it is easily seen that - compared to the zero- value 
minimum in the quantum regime - in the classical case R^ n \0) = 1/2 is valid. 

Let us now consider the case when a Bragg mirror is introduced in one arm of 
the interferometer (say, the second arm). In Fig. 3 the dependence on frequency 
of the square of the absolute value of the transmittance of an assumed barrier, 
|T 2 i(c<;)| 2 , and the absolute value of the band-shape function of an assumed 
time-limited wave packet, |/(u;)|, are plotted, |/(a>)| being centered at the 
middle of the band gap (2.68 PHz). Note that in the Berkeley experiments 
[8] the number of layers the barrier is composed of is N = 11. Examples of 
R( n \s) for the quantum regime are shown in Fig. 4 and discussed in detail in 
Ref. [10]. Here it should be pointed out that when N is not too large (e.g., N 
= 11 in Fig. 4) R^ is always smaller then unity, whereas for sufficiently large 
N (e.g., TV = 41 in Fig. 4) it becomes also possible that R^ > 1. Clearly, such a 
behaviour cannot be explained from interference of probability amplitudes of 
indistinguishable quanta, but this is not very surprising if we recall the strong 
distortion of a wave packet at such a barrier [10]. 

Since with increasing number of layers, N, the transmittance of the barrier 
drastically decreases, it may be difficult to observe the effect in the scheme 
considered so far. However, it may easily be observed if the laser is detuned 
such that the band-shape function f(uo) bridges over a band-gap edge of the 
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barrier and sufficiently overlaps a region of (fractional) transparency. From the 
interference fringes shown for the quantum regime in Fig. 5, in which |/(u;)| 
is assumed to be centered at the upper band-gap edge (3.06 PHz), it is clearly 
seen that the number of coincidence events can be reduced (R^ < 1) as well as 
enhanced (R^ > 1). Note that with decreasing N the main minimum survives, 
whose position at s < indicates subliminal behaviour. 

It could be thought that the effect of enhanced coincidences results from the 
distortion of the wave packet at the barrier relative to the undisturbed wave 
packet in the other arm of the interferometer. That this is not the case can be 
seen from Fig. 6 (for the quantum regime). In the figure it is assumed that the 
two wave packets are distorted in the same way, introducing identical Bragg 
mirrors in the two arms of the interferometer (see Fig. 1). Again, reduced and 
enhanced correlations are observed. In contrast to Fig. 5 however, R^(s) is 
now a symmetric function of s, R^(—s) = R^(s), because of the symmetry 
of the equipment. 

For comparison, let us now consider disentangled two-photon states in the 
quantum regime or two uncorrelated classical pulses, i.e., light beams which 
are produced by independent sources. In this case in Eqs. (9) and (12), re- 
spectively, ifi(s) [Eq. (11)] and Gi(s) [Eq. (14)] read as 



Ms) = 



J dQa(Q) J dujuj\f(uj)\ 2 e 2i<JJs/c T^\uj)T^\uj) 



(19) 



and 



G 1 ( 



' dujuje 2iu,s / c \f(uj)\' 




(20) 



It is seen that the number of coincidence events can only be reduced, i.e., 
R^(s) < 1, but not enhanced, R^(s) ^l. 

Enhanced coincidences (as shown in Figs. 4-6) can therefore be regarded as 
being a signature of correlated beams. However it also depends on the band- 
shape functions of the interfering wave packets, whether enhanced coincidences 
can be observed or not. To be more specific, let us consider the band-shape 
functions (u) = T 2 ( { /n) (u)f(u) of the wave packets in the two input ports 
of the beam splitter. The wave packets f^ u \u) will be said to be of type A 
if the Fourier transform 



oo 

f(4s) = J du [e Aisu)/c [e(w + \n)-e{cj- \n) 
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x (Itf - w ') f®( u ) (/Cn>( w ))* fW(-u>) (21) 



attains non-negative values only. Otherwise, the wave packets will be said to 
be of type B. In particular, it is easily seen that when |/(u;)| is a Gaussian and 
T 2 ? /n) (w) = l,then are of type A, whereas in the presence of a Bragg 

mirror, T^ 11 \uj) ^ 1, wave packets of type B can be produced. Recalling 
Eqs. (11) and (17), we find that for wave packets of type A the effect of 
enhanced coincidences cannot be observed on principle. 



4 Summary and conclusions 

We have studied the influence of wave-packet distortion at multilayer dielectric 
barriers on the coincidence-event output statistics in two-beam interference 
experiments of the type described in [8]. Usually the observed dependence 
of the interference fringes on the optical-path difference in the interferom- 
eter indicates that the wave packets prefer to arrive at equal output ports 
rather than different ones which in a sense can be interpreted as a "bunching" 
behaviour. However, when correlated wave packets are sufficiently distorted, 





Uncorrelated 
beams 


Correlated beams 


Band shape 


Type A or B 


Type A 


TypeB 


Classical 
beams 


R^\s) < 1 


i?( n )(s) < 1 
R^(s) > 1 
V<\ 


Single-photon 
quantum 
beams 


R^\s) < 1 
V < 1 


i?( n )(s) < 1 

R^\s) > 1 
V < 1 



Table 1 

Conditions of reduction and enhancement of the number of coincidence events in 
the two output ports of the beam splitter in the Mach-Zehnder interferometer in 
Fig. 1. 

then the opposite effect can also be observed, i.e., different output ports are 
preferably used and hence the number of coincidence events is enhanced, in 
correspondence to an "antibunching" behaviour of the wave packets at the 
beam splitter. The conditions under which the number of coincidence events 
is reduced or enhanced are summarized in Table 1, with special emphasis on 
the differences between quantum and classical light. 
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It should be noted out that similar effects are expected to be observed when 
the distortion of the wave packets is caused by other than multilayer devices. 
Finally, it should be pointed out that the diminished fringe visibility of the 
interference fringes in the classical regime has been derived for pulses produced 
in a random way in close similarity to the quantum description of the down- 
conversion process. If there is an appropriate correlation between consecutive 
wave packet pairs, this visibility may be enhanced up to 100% also in the 
classical situation [16]. 
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Figure 1: Scheme of an extended two-photon interference experiment of the 
type [8] for determining photon traversal times through multilayer dielectric 
barriers (L, laser; P, prism; DB I/II, dielectric barriers; BS, beam splitter; 
PDi, PD 2 , photodetectors). 
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Figure 2: The (normalized) coincidences R^ n \s) are shown for an experimental 
scheme as shown in Fig. 1 in dependence on the translation length s for time- 
limited pulses (2t = 40 fs) of photons traveling in free space. 
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|T2i(w)| , |f(^)| (normalized) 




uj/PHz 



N=ll N=57 ----- \f( u )\ 



Figure 3: The square of the absolute value of the transmittance of a mul- 
tilayer non-absorbing barrier (n Ti02 = 2.22, n Si02 = 1.41), |T 2 i(cj)| 2 , is shown 
for N — 11 layers and N = 57 layers together with the (normalized) spectral 
line shape function of a pulse that is assumed to be limited in time (2t = 
40 fs). 
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2s /c (fs) 



Figure 4: The (normalized) coincidences R^ n \s) are shown for an experimental 
scheme [8] in dependence on the translation length s for a time-limited pulse 
of the incoming photon (2t = 40 fs) and various numbers of the layers of a 
lossless barrier: TV = 11 (dotted-dashed line), iV = 35 (full line), TV = 41 (dashed 
line). The data of the lossless barrier are the same as in Fig. 3. 
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Figure 5: The (normalized) coincidences R^(s) are shown for an experimental 
scheme similar to [8] but for a laser frequency Q = 6.22 PHz, in dependence 
on the translation length s for a time-limited pulse of the incoming photon 
(2t = 40 fs) and a Bragg mirror with N = 57 layers. The data of the lossless 
barrier are the same as in Fig. 3. 
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Figure 6: The (normalized) coincidences R^(s) in dependence on the trans- 
lation length s are shown for the experimental scheme sketched in Fig. 1 
including two identical Bragg mirrors with N = 57 layers in path I and II. A 
laser frequency tuned to Q = 6.16 PHz and time-limited pulses of the incom- 
ing photons (2t = 40 fs) have been assumed. The data of the lossless barriers 
are the same as in Fig. 3. 
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